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Abstract 


The decay K~ —> n^e~u has been studied using in-flight decays detected with ’TSTRA-t-” 
setup operating at the 25 GeV negative secondary beam of the U-70 PS. About 130K events 
were used for the analysis. The A_|_ parameter of the vector formfactor has been measured: 
A+ = 0.0293±0.0015(stat)±0.002(s|/st). The limits on the possible tensor and scalar couplings 
have been derived: 

/r//-i-(0) = —0.044lgQ57 (stat) 
fs/hW = -0.020in?i (stat) 
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Figure 1: The layout of the ’TSTRA+” setup. 

1 Introduction 

The decay K —»• evTT^{Ke^) is known to be a promising one to search for an admixture of scalar 
(S) or tensor (T) interactions to the Standard Model (SM) V-A . This topic has been attracting 
significant interest during resent years and moreover some previous experiments with charged 
and neutral kaon beams have reported indications for some anomalous S and T signals |l| , [@] . 
On the other hand, resent KEK [Q, Q experiment with stopped beam has reported negative 
results of the searches. 

In our analysis we present new search for S and T couplings based on the statistics of about 
130K Ke3 events. Another result of our study is the measurement of the V-A /+(f) formfactor 
slope A+. 

2 Experimental setup 

The experiment is performed at the IHEP 70 GeV proton synchrotron U-70. The experimental 
apparatus ’TSTRA-I-” is the result of the modification of’TSTRA-M” |^, which, in turn, evolved 
from ’’ISTRA” that yielded important results on ti~ and K~ decays in the late 1980’s |^. The 
setup is located at the 4A negative unseparated secondary beam. The beam momentum in the 
measurements is ~ 25 GeV with Ap/p ~ 2%. The admixture of K~ in the beam is ~ 3%. The 
beam intensity is ~ 3 ■ 10® per 1.9 sec. U-70 spill. A schematic view of the detector is shown 
in Fig.l. The momentum of the beam particle deflected by Mi is measured by BPCi BPC^^ 
PG’s with 1mm wire step, the kaon identihcation is done by Ci-^C^ threshold G-counters. The 
9 meter long vacuumated decay volume is surrounded by 8 lead glass rings LGi -=- LG^ used 
to veto low energy photons. The same role is played by SP 2 - a 72-cell lead glass calorimeter. 
The decay products deflected in M2 with ITm field integral are measured with PGi -i- PG^- 
2mm step proportional chambers; DGi -h DG^- 1cm cell drift chambers and finally with 2cm 
diameter drift tubes DTi -h DT^. A wide aperture threshold Gerenkov counter C 4 is filled 
with He and used to trigger the electrons. SPi is a 576-cell lead glass calorimeter, followed by 
HG - a scintillator-iron sampling hadron calorimeter, subdivided into 7 longitudinal sections 
7x7 cells each. MH is a 11x11 cell scintillating hodoscope, used to solve the ambiguity for 
multitrack events and improve the time resolution of the tracking system, MuH is a 7x7 cell 

















































muon hodoscope. 

The trigger is provided by Si S '5 scintillation counters, Ci C 3 Cerenkov counters, analog 
sum of amplitudes from the last dinodes of the SPi and is very loose: T = Si ■ S 2 ■ S 3 ■ S 4 ■ Ci ■ 
C 2 ■ C 3 ■ ^5 ■ S(S'Pi), here S4 is a scintillator counter with a hole to suppress beam halo ; S'5 is a 
counter downstream the setup at the beam focus; S(SPi) - a requirement for the analog sum of 
amplitudes from SPi to be larger than ~700 MeV - a MIP signal. The last requirement surveys 
to suppress the iC —>• /iz/ decay. Some complementary triggers: Tg = S'! ■ S '2 ■ S '3 • S '4 • S '5 ■ (T 4 - the 
electron trigger and prescaled ” decay” trigger Td = Si ■ S 2 ■ S 3 ■ S 4 ■ S 5 were used to cross-check 
the efficiency of the main one. 

The main difference between ’TSTRA-M” and ’TSTRA-I-” is in the electronics and DAQ: 
all the CAMAC based electronics was changed by IHEP developed MICC [0 ECL-based elec¬ 
tronics. ’’ISTRA-I-” has now 12 MICC crates with ADC’s, TDC’s and latches. The DAQ, 
described in some details in [§] is based on IHEP-developed VME master V-08 0 , which writes 
the MICC stream into standard VME memory. Between the spills, the information is written 
into PC through BIT-3 VME-PCI interface. The saturated event rate is ~ 6500 of 1 Kb events 
per 1.9 sec. spill. 


3 Event selection 


During physics runs in November-December 1999 and March-April 2001, 206M and 363M 
events were logged on DLT’s. This information is supported by about lOOM MC events gener¬ 
ated with Geant3 |^. MC generation includes realistic description of the setup: decay volume 
entrance windows, track chambers windows, gas, sense wires and cathode structure, Cerenkov 
counters mirrors and gas, shower generation in EM calorimeters etc. 

The usual first step of the data processing is the EM calorimeter calibration, using special 
runs with 10 GeV electrons; track system alignment, HCAL and guard system calibration 
with muon beam runs. The data processing starts with the beam particle reconstruction in 
BPCi-rBPC/i, then the secondary tracks are looked for in PC'i-i-PC '3 ; PCi-i-PCs; DTi^DT 4 
and events with one good negative track are selected. The decay vertex is searched for, and 
a cut is introduced on the matching of incoming and decay track. The next step is to look 
for showers in SP 2 calorimeter. A method of shower parameters reconstruction based on the 
MC-generated patterns(~ 2000 3x3 patterns) of showers is used. The matching of the charged 
track and a shower in SP 2 is done on the basis of the difference between the track extrapolation 
and the shower coordinates. The electron identification is done using E/P ratio - of the energy 
of the shower associated with the track and the track momentum, see Fig.2. The selection 
of events with the two extra showers results in spectrum shown in Fig.3. The vr® peak 
has a mass of M^ro = 134.8MeV, and a resolution of 8.6 MeV. Another important variable 
for the K~ —»• ez/ 7 r° selection is the missing mass squared: {Pk — Pe — P-k^Y-, where P are the 
corresponding four-momenta, see Fig.4. The cut is ±0.01 GeV^. The further selection is done 
by the requirement that the event passes 2 C K —y fit. At the same time, similar 2C 
fit A' —>■ TT~Ti^ should fail. The missing energy — E^o after this selection is shown 

in Fig.5 The peak at low E^iss corresponds to the remaining K~ —>■ 7 r“ 7 r° background. The 
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Figure 2: The E/p plot - the ratio of the energy of the associated cluster in ECAL to the 
momentum of the charged track. The arrow shows the cut used for the electron separation. 
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Figure 3: The 77 mass spectrum for the events with the identified electron and two extra 
showers. 
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Figure 4: The missing four-momentum squared {Pk — Pe — P-k^Y for the selected events for 
1999 run (left) and 2001 run (right). The points with errors are the data, the histogram - MC. 

corresponding cut is Emiss > ^GeV. The surviving background is estimated from MC to be 
less than 3%. 

The detailed data reduction information is shown in Table. 1. 


4 Analysis 


The event selection described in the previous section results in selected 54K events in 1999 
data and 79K events in 2001 data. The distribution of the events over the Dalitz plot is shown 
in Fig.6. The variables y = 2 E(./Mk and 2 ; = 2E-,^/Mk-, where E^^ Et^ are the energies of the 
electron and 7r° in the kaon c.m.s are used. The background events, as MC shows, occupy the 
peripheral part of the plot. 


The most general Lorentz invariant form of the matrix element for the decay K —> I 


IS 
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M = 2f:^S(pJ(l + 7")|m/f/s - \[(Pk + + (Pk - + t—a^^P‘PS]v{pi) 

V 2 ^ 'mR 

( 1 ) 

It consists of scalar, vector and tensor terms, fs, fr, f± are the functions of f = {Pk — P-kY- In 
the Standard Model (SM) the W-boson exchange leads to the pure vector term. The ’’induced” 
scalar and/or tensor terms, due to EW radiative corrections are negligible small, i.e the nonzero 
scalar/tensor form factors indicate a physics beyond SM. 

The term in the vector part, proportional to /_ is reduced (using Dirac equation) to a scalar 
formfactor. In the same way, the tensor term is reduced to a mixture of a scalar and a vector 
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Figure 5: The missing energy for the evr^ events, for 1999 run (left) and 2001 run (right). 
The points with errors are the data, the histograms - MC. The dark(blue) peak at zero value 
corresponds to the MC-predicted K —>■ background. The arrow indicates the cut value. 
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Figure 6: Dalitz plots {y = 2Ee/MK]z = 2 Et,o/Mk) for the selected K —>• ez/7r° events after 
the 2-C £t. Left - 1999 statistics. Right - 2001 statistics. 
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Table 1: Event reduction statistics. The main steps are shown for the 1999 and 2001 runs . 


Run 

1999 

2001 

Nevents OU tapeS 

206.544.909 

363.002.105 

Beam track reconstructed 

159.459.629=77% 

268.564.958 =74 % 

One secondary track found 

81.166.929=41% 

134.227.095 =37% 

Written to DST 

70.015.610 =34% 

107.215.783 =30 % 

e~ identihed 

1.300.958 

1.998.719 

2 showers reconstructed 

252.177 

361.621 

7 r° identihed 

186.850 

251.489 

WLJ < 0.01 

96.652 

144.642 

K —>■ 7i~71^ rejected 

79.660 

117.566 

K —>• accepted 

65.208 

97.585 

Emiss > 1 GeV 

54.009 

79.248 


formfactors. The redehned /+(V), ^^(S) and the corresponding Dalitz plot density in the kaon 
rest frame(p(E 7 r, E;)) are [|I2| : 


V = 

S = 

p(B^, El) ~ 
A = 

B = 
C = 
= 


f+ + {rni/mK)fT 

fs + {mi/2mK)f- + - 

V 

A - \V\^ + B-Re{V*S) + C -ISl^ 

TTiKi^EiEy - ttikAE^) - m^i{E^ - 
mimK{2E^ - AE^) 

^max _ ^ 

TT 


2Ee _ E^ \ 
rriK rriK) 



/r 


( 2 ) 


In case of Ke3 decay one can neglect the terms proportional to mi] mf. Then, assuming linear 
dependence of /+ on t: /+(t) = /+(0)(1 + and real constants fs, fr we get: 


p{E^, El) 


rsj 


+ 


mK{2EiE^ 
m\AE^ ■ ( 


- mxAE^) ■ (1 + X+t/mlf 

fs , [i_ 2 -^e _ /t 

/+( 0 ) V mx mKJ /+( 0 ) 


2 


(3) 
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The procedure for the experimental extraction of the parameters A+, fs, /t starts from the 
subtraction of the MC estimated background from the Dalitz plots of Fig.6. The background 
normalization was determined by the ratio of the real and generated K~ 7r“7r° events. Then 
the Dalitz plots were subdivided into 20 x 20 cells. The background subtracted distribution of 
the numbers of events in the cells (i,j) over Dalitz plots, for example, in the case of simultaneous 
extraction of A+ and was fitted with the function: 

p(b J) ~ + W2it,j) ■ A+ + ■ A^ + W,ii,j) ■ (4) 


Here Wi are MC-generated functions, which are build up as follows: the MC events are generated 
with constant density over the Dalitz plot and reconstructed with the same program as for the 
real events. Each event carries the weight w determined by the corresponding term in formula 
3, calculated using the MC-generated values for y and z. The radiative corrections according to 


T3[| were taken into account. Then Wi is constructed by summing up the weights of the events 
in the corresponding Dalitz plot cell. This procedure allows to avoid the systematics errors 
due to the ’’migration” of the events over the Dalitz plot because of the finite experimental 
resolution. 


5 Results 


The results of the fit are summarized in Table.2. The combination of the two runs is done by 
the simultaneous fit. The first line corresponds to pure V-A SM fit. In the second line the 
tensor and in the third the scalar terms are added into the fit. All the errors presented are 
from the ’’MINOS” procedure of the ’’MINUIT” program [M| and are larger than the Gaussian 
ones. At present, we estimate an additional systematics error in A+ to be ±0.002. The estimate 
is done by comparing two runs, which differ a lot in amount of matter in the beamline and 
detector conhguration and by varying cuts, cell size during the £t of the Dalitz plots etc. 

The comparison of our results with the most recent data [2, ^ shows that our 

statistics, at present, is the highest in the world and the statistical errors somewhat smaller 
than in []^, and comparable with [Q. We do not confirm the observation of a significant fs 
and /t in ||I|]. Our data is in a good agreement with [§, ^ and with the theoretical calculations 
for A_|_: A_|_ = 0.031 [^, done in the context of the chiral perturbation theory. 


6 Summary and conclusions 

The decay has been studied using in-flight decays of 25 GeV K~, detected by ”ISTRA±” 
magnetic spectrometer. Due to the high statistics, adequate resolution of the detector and 
good sensitivity over all the Dalitz plot space, the measurement errors are significantly reduced 
as compared with the previous measurements. The A+ parameter of the vector formfactor has 
been measured to be: 


A+ = 0.0293 ± 0.0015(stot) ± 0.002(s|/st). 


















Table 2: Results of the fit. 



1999 

2001 

1999+2001 

A+ 

n no 7 i +0.0023 
U.UZ ( 1_0.0023 

0 031 

U.UOiU_o.ooi9 

0.02931°:°°}® 


A+ 

n 0970 + 0.0023 
U.UZ /U_o,o023 

0.0310+“19 

0.02931°:°°}® 

/t 

= Fr/f+m 

1 

0 

0 

00 

00 

00 

1 + 

0 0 

b b 

00 00 

00 00 

-0.04871°:°?°® 

-0.0445l°:°f7l 


A+ 

0.0268;°;“27 

0 n304'''°'°°^^ 

U.UOU+_o.o024 

0.02891°:°°}^ 

fs 

= Fs/UiO) 

-0.0139+°ii 

0 O99c:+ 0.0355 
—U.UZZO_o.oi87 

-0.0i97l°:°?i 

xVndf 

1.7 

1.3 

1.5 

"A^bins 

225 

228 



The limits on the possible tensor and scalar couplings have been derived: 

/t//+(0) = -0.0441^^?; 

fs/f+{0) = -0.020tH?i 
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